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Abstract 

Radiatively generated, LO quark (u, d, s, c, b) and gluon densities in the real, unpo- 
larized photon, calculated in the CJK model being an improved realization of the CJKL 
approach, have been recently presented. The results were obtained through a global fit to 
the experimental i^2^(a;, Q^) data. In this paper we present, obtained for the very first time 
in the photon case, an estimate of the uncertainties of the CJK parton distributions due to 
the experimental errors. The analysis is based on the Hessian method which was recently 
applied in the proton parton structure analysis. Sets of test parametrizations are given 
for the CJK model. They allow for calculation of its best fit parton distributions along 
with F2{x,Q'^) and for computation of uncertainties of any physical value depending on 
the real photon parton densities. We test the applicability of the approach by comparing 
uncertainties of example cross-sections calculated in the Hessian and Lagrange methods. 
Moreover, we present a detailed analysis of the of the CJK fit and its relation to the 
data. We show that large x^/dof of the fit is due to only a few of the experimental 
measurements. By excluding them x^/dof ~ 1 can be obtained. 
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1 Introduction 



In the recent paper, [Ij, new results on the LO unpolarized real photon parton distributions 
have been presented. In that work we improved and broadened our previous analysis, j2]. 
As a result three models and corresponding parametrizations of the photon structure function 
F2{x,Q'^) and parton densities were given. 

The main difference between the models lies in the way the heavy charm- and bottom-quark 
contributions to the photon structure function are described. In two approaches, referred to as 
FFNScjii-1 and 2 we adopted a widely used scheme in which there are no heavy quarks, h, in 
the photon. In that case heavy quarks contribute to the photon structure function through the 
so-called Bethe-Heitler, 7*7 — hh interaction. The FFNScja'2 model includes an additional 
"resolved" contribution given by the 7*^ —>■ hh process, see for instance P]. Finally, the model 
of our main interest, denoted as CJK, applies the ACOT(x) |1] scheme, where heavy-quark 
densities appear. 

The free parameters of each model have been computed by the means of the global fits to the 
set of updated data collected in various e~^e~ experiments. That way parametriza- 

tions of the photon structure function and parton distributions were created. 

The main goal of the present analysis is to estimate the uncertainties of the CJK parton 
distributions due to the experimental errors of F2{x, Q"^) data. Alike in the global CJK fit the 
raw experimental data are used, by which we mean that neither the radiative corrections, nor 
the corrections taking into the account the small off-shelness of the probed quasi-real photons 
are included. Further, the correlations among the measurements are neglected and in each x bin 
the approximation ~ is assumed. Finally, one has to keep in mind that all the 

data were obtained with an assumption of the FFNS scheme for the heavy-quark contributions 
to -F^(a;, Q"^) . The inclusion of the above corrections is beyond the scope of this work. 

Our work has been motivated by the recent analysis performed for the proton structure 
by the CTEQ Collaboration, |S]-|Zj and the MRST group, jH]. We use the Hessian method, 
formulated in recent papers, to obtain sets of parton densities allowing along with the parton 
distributions of the best fit to calculate the best estimate and uncertainty of any observable 
depending on the photon structure. As the same aim can be obtained utilizing the Lagrange 
method, see for instance [Sj,[Sl and P, we apply this approach to obtain an independent test 
of the correctness of the Hessian model results. 

Moreover we perform an analysis of the sources of the obtained in the global fit for CJK 
model. Namely, we examine which e~^e~ experiments data produce largest contributions to the 
total of this fit as well as for other models including our FFNScj_ft:l, 2 approaches and GRS 
LO and SaSlD fTT| parametrizations. 

The paper is divided into five parts. Section 2 shortly recalls the FFNScji<- and CJK 
models of the real photon structure. Section 3 is devoted to an introduction of the Hessian 
method and to a presentation of the calculated uncertainties of the CJK parton distributions. 
Furthermore, in section 4, we briefiy introduce the Lagrange method. Next, in section 5, 
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we compare the uncertainties of two example cross- sect ions calculated in both, Hessian and 
Lagrange approaches. Finally, in section 6 we examine the data sources of the of the CJK 
fit. The parton distributions which are a result of our analysis have been parametrized on the 
grid. The FORTRAN programs obtained that way are open to the public use. 

2 FFNS(:7jir and CJK models - short recollection 

In this analysis we focus on the CJK model results. We are interested in the parton distributions 
computed in this approach and most of all in their uncertainties due to the experimental data 
errors. Still, in order to estimate the allowed deviation of the of the global CJK fit from its 
minimum we will apply the FFNS(7j;^ models. Moreover, they will be very useful in determining 
of the e"'"e~-experiment data sets which are not well described by the photon structure models. 
Therefore, we shortly recall the main differences between the FFNS type approaches and the 
CJK approach. Next, we present the origin of the free parameters of the CJK model calculated 
in jlj by the means of the global fit to the -^^(x, Q"^) data. Further, the data used in the fits is 
described. Finally, we recall the results of the CJK model fit. 

2.1 The models 

The two approaches leading to the FFNSpji^- and CJK models, have been described in detail 
in our previous papers PP and [2]. The difference between them lays in the approach to the 
calculation of the heavy, charm- and beauty-quark contributions to the photon structure func- 
tion Q^). First, FFNScji^- type models base on a widely adopted Fixed Flavour Number 
Scheme in which there are no heavy quarks in the photon. Their contributions to F2{x,Q'^) 
are given by the 'direct' (Bethe-Heitler) 7*7 —>■ hh process of heavy-quarks production h. In 
addition one can also include the so-called 'resolved'-photon contribution: 7*6* ^ hh. We 
denote these terms as F2i^{x,Q^)\direct and i^{x, Q'^)\resoived, respectively. We considered two 
FFNS models: in the first one, FFNScji^l, we neglected the resolved photon contribution, 
while in the second one, FFNS(7jft'2, both mentioned contributions to Q^) were included. 
In such models the heavy-quark masses are kept to their physical values. The photon structure 
function is then computed as 

3 

F2{x,Q')=J2^eK^l + ^l)i^^Q')+ E [FU^,Q')Urect + Fl,{x,Q')\resol.e,] , (1) 
i=l h{=c,b) 

with qJi^x^Q"^) {q]{x,Q^)) being the light u,d,s quark (anti-quark) densities, governed by the 
Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) evolution equations. 

The CJK model adopts the new ACOT(x) scheme, 0, which is a recent realization of the 
Variable Flavour Number Scheme (VFNS). In this scheme one combines the Zero Mass Variable 
Flavour Number Scheme (ZVFNS), where the heavy quarks are considered as massless partons 
of the photon, with the FFNS just discussed above. In this model, in addition to the terms 
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shown in Eq. (P) one must include the contributions due to the heavy-quark densities which 
now appear also in the DGLAP evolution equations. A double counting of the heavy-quark 
contributions to F2{x,Q'^) must be corrected with the introduction of subtraction terms for 
both, the direct and resolved-photon, contributions. Further, following the ACOT(x) scheme, 
we introduced the Xh = x{l + 4m1/Q'^) parameters responsible for the proper vanishing of the 
heavy-quark densities at the kinematic thresholds for their production in DIS: 
W'^ = Q'^{l—x)/x > 4:m\, where W is the 7*7 centre of mass energy. Adequate substitution of x 
with Xh in Qh and the subtraction terms forces their correct threshold behavior as Xfc 1 when 
W — > 2mh- This is achieved for all the terms except for the subtraction term Q^)\dir,subtr 
which requires an additional constraint, F2 f^{x,Q'^)\dir,subtr = for Xh > 1, imposed by hand. 
Finally, we obtain the following formula for the photon structure function in the CJK model 

F^ix,Q') = xj^el{ql + ql){x,Q') + xY.el{ql + ql){xH,Q') 

i=l h{=c,b) 

~l~ ^ ^ ['^2./i ('^' ^5 ) Idfrerf ~l~ -^2,/i("''' Q ) | resoived] (2) 
h{=c,b) 

h(=c,b) 



The above formula for Q^) must be complited by imposing of another condition (the 

positivity constraint): 

F2,h{-^iQ ) — -^2^/1, (-^i Q ) I direct + -^2^/1(3;, Q )\resolved (3) 

preventing the unphysical situation F^f^{x,Q'^) < F^^^{x,Q'^)\direct + F^fX^,Q'^)\resoived which 
can occur at small- and large-x. 

Explicit expressions for the terms appearing in Eqs. and can be found in [2j. 

For all models we chose to start the DGLAP evolution at a small value of the scale, 
Qo = 0.25 GeV^, hence our parton densities are radiatively generated. As it is well known the 
point-like contributions are calculable without further assumptions, while the hadronic parts 
need the input distributions. For that purpose we utilize the Vector Meson Dominance (VMD) 
model ^2); with the assumption that all light vector meson contributions are proportional 
to the p° meson contribution and are accounted for via a parameter k, that is left as a free 
parameter. We took parton densities in the photon equal to 

rHaMl) = ^^rMl)- (4) 

The parameter is extracted from the experimental data on the r(p° — > e+e^) width. 

We assumed the input densities of the p° meson at Ql = 0.25 GeV^ in the form of valence-like 
distributions both for the (hght) quark (v^) and gluon (C) densities. All sea-quark distributions 
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(denoted as (p), including s-quarks, were neglected at the input scale: 



xGP{x, Ql) = Ngxv^ix, Ql) = Ngx^'il - x)^, 
xC{x,Ql) = 0, 



(5) 
(6) 
(7) 



where A^^ = NgN^. 

The valence-quark and gluon densities must satisfy the constraint representing the energy- 
momentum sum rule for p: 



Jo 

That constraint allowed to express in the case of the CJK model the normalization factor A^^^ 
as a function of a, f5, N^, and k. 

In this analysis we decided to relax the second possible constraint, related to the number 
of valence quarks in the p° meson, n^: 



We observed that imposing that constraint spoils the quadratic approximation of the Hessian 
method described in section 01 The reason may be the additional correlation between the a 
and (3 parameters which appears when we express the Ny parameter as a function of the above 
two. On the other hand, in the CJK model case, the global fit without imposing the additional 
constraint gives proper value, = 2 ± 0.1. Therefore we decided to use it as the CJK fit 
and further apply the Hessian method to examine the uncertainties of the resulting parton 
distributions. 

2.2 Data 

The CJK as well as the FFNScjir model fits were performed using all existing Q^) 
data, |ini^|21I, apart from the old TPC27, [22] • In our former global analysis [2j we used 
208 F^(a;, Q^) experimental points. Now we decided to exclude the TPC27 data from the set 
because it has been pointed out that these data are not in agreement with other measurements 
(see for instance [IE])- After the exclusion we are left with 182 F2{x,Q'^) experimental points. 

The raw experimental data are used, by which we mean that neither the radiative cor- 
rections, nor the corrections taking into the account the small off-shelness of the quasi-real 
photons probed in the experiments are included. Further, the correlations among the measure- 
ments, not presented in most of the articles 5ni^[21], are neglected. Finally, in each x bin the 
approximation (F2^(x)) ~ F2{{x)) is assumed. 

We included all the data in the fit without any weights. A list of all experimental points 
used can be found on the web-page fT7\ . 
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2.3 Results of the CJK global fit 



The fit of tlie CJK model to tlie experimental F2{x,Q'^) data set described above was based 
on tlie least-squares principle (minimum of x^) and were done using MiNUlT ^2^. Systematic 
and statistical errors on data points were added in quadrature. 

The results of our new fit are presented in table ^ The second and third columns show the 
quality of the fit, i.e. the total ^oi 182 points and the P^r degree of freedom. The fitted 
values for parameters a, (3, k and A^^, are presented in the middle of the table with the errors 
obtained from MiNOS with the standard requirement of Ax^ = 1. In addition, the value for Ng 
obtained from other parameters using the constraint (jH)) is given in the last column. 



model 


(182 pts) 


X^/dof 


K 


a 








CJK 


273.7 


1.537 


1 OQA+O.lSl 


U.Z»y_Q Qgg 


0.898+°i^^ 


U.^U^_Q Qgg 


4.93 



Table 1: The total for 182 data points used in the fit and for the degree of freedom and 
parameters of the fit for the CJK model. All the given errors are obtained from MiNOS with 
the standard requirement of Ax^ = 1. 



The results of the CJK fit, its agreement with various experimental data and its comparison 
to the other models and other photon structure parametrizations have been described in detail 
in . In this paper we focus on the uncertainties of the parton distributions due to the exper- 
imental F2{x, Q^) data. Moreover, we check in detail the influence of various measurements on 
the quality of the fit. 

3 Uncertainties of the parton distributions - Hessian 
method 

We are going to address in this article a problem that so far has never been considered in 
the case of the photon structure. Namely, we want to present an analysis of the experimental 
uncertainties of our CJK parton densities. To reach that goal we derive the necessary knowledge 
and tools from the proton structure considerations. 

During the last two years numerous analysis of the uncertainties of the proton parton densi- 
ties resulting from the experimental data errors appeared. The CTEQ Collaboration in a series 
of publications, I5]-[7j, developed and applied a new method of their treatment significantly 
improving the traditional approach to this matter. Later the same formalism has been applied 
by the MRST group in |H]. The method bases on the Hessian formalism and as a result one 
obtains a set of parametrizations allowing for the calculation of the uncertainty of any physical 
observable depending on the parton densities. 
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The Hessian method is a very useful tool as it allows for computing of the parton density 
uncertainties in a very simple and effective way. Still, it relies on the assumption of the quadratic 
approximation which not necessary is perfectly preserved. 



3.1 The method 

A detailed description of the method applied in our analysis can be found in |5j and jj^. For the 
sake of clearness of our procedures we will partly repeat it here, keeping the notation introduced 
by the CTEQ Collaboration. 

Let us consider a global fit to the experimental data based on the least-squares principle 
performed in a model, being parametrized with a set of {a^, i = 1,2, ■ ■ ■ ,d} parameters. Each 
set of values of these parameters constitutes a test parametrization 5*. The set of the best values 
of parameters corresponding to the minimal x^? Xo' denoted as parametrization. In 
the Hessian method one makes a basic assumption that the deviation of the global fit from Xo 
can be approximated in its proximity by a quadratic expansion in the basis of parameters {aj} 

d d 

Ax' = - Xo = E E ^^^("^ - «°)(«^- - (10) 
i=i j=i 

where Hij is an element of the Hessian matrix calculated as 



• ^^^^ 



ao 



Since the Hij is a symmetric matrix, it has a complete set of k = 1, 2, ■ ■ ■ ,d orthonormal 
eigenvectors {vi)k defined by 



^Hij{vj)k = tk{vi)k, (12) 
i=i 

d 

^iyi)j{vi)k = Sjk, (13) 



i=l 



with {efc} being the corresponding eigenvalues. Variations around the minimum can be ex- 
pressed in terms of the basis provided by the set of eigenvectors 



ttj — a° 



'^SkZk{vi)k, (14) 



fc=i 



where {zk} are new parameters describing the displacement from the best fit. The {sk} are 
scale factors introduced to normalize {zk} in such a way that 

AX' = X'-X^ = X:4 (15) 

k=l 
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The above equation means that the surfaces of constant are spheres in {z^} space. That way 
the {^z\~\ coordinates create a very useful, tenormahzed basis. The {vi)k = matrix describes 
the transformation between this new basis {zk} and the old {oj} parameter basis. The scaling 
factors Sk are equal to ^/TJek provided that we work in the ideal quadratic approximation. In 
reality they differ from these values. 

The Hessian matrix can be calculated from its definition in Eq. (fTT|) . Such a computation 
meets many practical problems arising from the large range spanned by the eigenvalues {6^}, 
the numerical noise and non-quadratic contributions to x^- The solution has been given by 
the CTEQ Collaboration [Sj which introduced an iterative procedure working properly in the 
presence of all the problems listed above. The algorithm has been implemented as an extension 
to the MiNUlT program and is open for public use (see 0). 

Having calculated the eigenvectors, eigenvalues and scaling factors we can create a basis 
of the parametrizations of the parton densities, {S^,k = 1, ■ ■ ■ ,d}. Each of the pairs 
corresponds to a different eigenvector direction {vi)k- These parametrizations are constructed 
in the following way: their parameters are defined by displacements of a magnitude t "up" or 
"down" along the corresponding eigenvector direction 



For example the set of parameters is given by inserting {zi, Z2, z^,- ■ ■) = (t, 0, 0, ■ ■ ■ ) into 
Eq. (|14p. For each parametrization Ax^ = t^- 

Further, let us consider a physical observable X depending on the photon parton distribu- 
tions. Its best value is given as X{S^). The uncertainty of X, for a displacement from the 
parton densities minimum by A^^ = (T - the tolerance parameter) can be calculated with 
a very simple expression (named as master equation by the CTEQ Collaboration) 



Note that having calculated AX for one value of the tolerance parameter T we can obtain 
the uncertainty of X for any other T by simple scaling of AX. This way sets of {S^} parton 
densities give us a perfect tool for studying of the uncertainties of other physical quantities. 
One of such quantities can be the parton densities themselves. 

Finally, we can calculate the uncertainties of the parameters of the model. According 
to Eq. (fT^ in this case ai(5'^) — cii{S^) = 2t{vi)kSk and the master equation gives a simple 
expression 



We end this shortened description of the used method with one practical point. In real 
analysis we observe the considerable deviations from the ideal quadratic approximation of 




(16) 




(17) 
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equation (jl5|) . To make an improvement we can adjust the scaling factors {sk} either to obtain 
exactly Ax^ = aX Zk = t for each of the 5*^ sets or to get the best average agreement over 
some Zk range (for instance for Zk <t). We chose to apply the second approach. 

3.2 Estimate of the tolerance parameter T for the CJK photon den- 
sities 

We consider now the value of the tolerance parameter T defining the allowed deviation of the 
global fit from the minimum, Ax^ = T^, as described in the previous section. Through the 
master equation (fTHj) T is relevant to the calculations of the uncertainties due to the real photon 
parton densities. In case of an ideal analysis Ax^ = 1 is a standard requirement. Of course 
a global fit to the F^(a;, Q^) data coming from various experiments is not such a case and 
certainly T must be greater than 1. Unfortunately no strict rules allowing for estimation of the 
tolerance parameter exist. A detail analysis of that problem can be found in [H] and [H]. We 
try to estimate the reasonable practical T value for the CJK fit in two ways, sa 

First we examine the mutual compatibility of the experiments used in the fit. We divide 
the data into six sets. Four of them contain the results of the CERN-LEP accelerator exper- 
iments - ALEPH, DELPHI, L3 and OPAL. The other two are sets of data collected with the 
DESY-PETRA and KEK-TRISTAN accelerators. The DESY-PETRA collection, referred to 
as DESY, combines the results of PLUTO, JADE, CELLO and TASSO collaborations. The 
KEK-TRISTAN set, denoted as KEK, combines the results of TOPAZ and AMY experiments. 
For each of the data sets we calculate the Xn and X-n ^^e nth collection), values of the of 
the best fit corresponding to the given experiment and to the remaining five ones respectively. 
Further we test how much x'^n can be lowered by minimazing with the removed nth set of 
data. We obtain Ax^„ value which is the minimal deviation of the global fit from its minimum 
necessary to describe the inclusion of the nth experiment to the global set of data. Results of 
that test are presented in table |21 The Xn/-^n values, where Nn is a number of the experimental 
points in the n-th data set, indicate that truly our global fit is not a case of an ideal analysis. 
For some of experiments our CJK fit agrees very well with the data. In case of others x^/Nn 
is much larger than 1 and x^/dof = 1-537 of the global CJK fit presented in tabled Further 
we see that the Ax^_^ varies from 0.1 to 20.9 obtaining the maximal value in the DELPHI 
experiments case. For full reliability of the test we performed an additional fit including all 
available -^^(x, Q^) data which means that apart from all the experiments used in the CJK 
fit and mentioned above the TPC27 data was utilized. We checked then that the Ax^_^ for 
the TPC27 case equals 3.2. We notice that the Ax^„ values in the case of DESY, DELPHI 
and OPAL sets are larger. One of the reasons for that may be the following: the Q^) 
points measured by the TPC27 experiment lie mostly in the low-Q^ region not covered by other 
measurements. That is not the case for the other experiments and so if their results differ the 
exclusion of one of them can strongly effect the x^ of the fit. Anyway, though there may exist 
other arguments for that, it seems difficult to support the statement that the TPC27 data are 
inconsistent with other measurements as claimed in [211 • Still, please notice that as described 
in section 12.21 the raw experimental data were applied and their various possible corrections 



9 



as well as the correlations among the measurements were not taken into consideration in our 
analysis. More discussion on the data will be given in the section IHl Finally, we have to assume 
that the allowed Ax^ is greater than 20.9 and hence the tolerance parameter must have value 
T ~ 5. 



n 


Set 


# of points Nn 


An 


XliNn 


A — n 




1 


DESY 


38 


89.5 


2.36 


184.2 


9.6 


2 


KEK 


16 


18.2 


1.14 


255.5 


0.2 


3 


ALEPH 


20 


21.0 


1.05 


252.7 


1.3 


4 


DELPHI 


38 


88.6 


2.33 


185.1 


20.9 


5 


L3 


28 


14.9 


0.53 


258.8 


0.6 


6 


OPAL 


42 


41.3 


0.98 


232.4 


5.8 



Table 2: The CJK model. Table presents number of the points, A^n in each data set, Xn (and 
Xn/Nn) and X-n^ being values of the x^ of the best fit corresponding to the nth set and to the 
remaining five sets respectively. Finally Axn is the value by which x'^n can be improved by the 
additional minimization. 

As a second test we compare the results of our three fits presented in this paper. We find 
the T values for which parton densities predicted by the FFNS models lie between the lines of 
uncertainties of the CJK model parton distributions. We made such test independently for each 
of the five quark flavours and for the gluon densities in three ranges and for 10^^ < x < 0.97. 
As can be seen in table IHl the resulting numbers differ for various parton densities. They also 
depend on the range in which we studied the T parameter. While not strong in the case 
of quark distributions this dependence is dramatic in the gluon case predicting large T value 
necessary to contain other gluon distributions within the CJK uncertainty-bands at high Q^. 

We do not think that the comparison of the Q^) values predicted by various models 

and parametrizations should give us better information on the necessary T parameter than the 
above test done for the parton distributions case. Q^) predictions seem to be much more 

model dependent than the resulting parton distributions. As an example we can mention the 
difference between the FFNScji^ fits (Ax^ ^ 34 - see table ^ which is not so apparent in 
the differences between the corresponding parton distributions. That can be easily understood 
when we recall that the FFNSc ji^- 1 and 2 models differ only by the resolved 'y*G ^ hh term 
which contributes to the Q^) but not to the DGLAP equations and therefore directly 

does not influence the parton densities. 

We estimate that the tolerance parameter T should lie in the range 5 ~ 10. Still, one has to 
keep in mind that because of the lack of data constraining the real photon gluon distribution, in 
case of processes dominated by the gluon interactions such assumption may not be safe enough. 
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T(G^) 




T(mT) 


T(sT) 




1 < < 100 GeV^ 


4.5 


7.0 


7.0 


3.4 


8 


1 < g2 < 1000 GeV2 


14.0 


7.0 


7.0 


3.4 


10.5 


1<Q'' < 200000 GeV^ 


138.0 


7.0 


7.0 


3.4 


20.0 



Table 3: The T[q^) values for which parton q'^{x,Q'^) densities predicted by the FFNS models 
lie between the lines of uncertainties of the CJK model parton distributions. Calculation 
performed for various ranges and for 10~^ < a: < 0.97. 

3.3 Tests of quadratic approximation 

As a result of application of the Hessian method to the CJK model we obtained a set of 
the {{vi)k} and {sk} values with i and k = 1, ■ ■ - 4 which corresponds to the number of free 
parameters (k, a, (3 and iV„) in the model. We used the iteration procedure provided by the 
CTEQ Collaboration (see The procedure was run with displacements giving Ax^ = 1, 3, 5 
and 10. Each time 15 iterations were performed. Tests of the quadratic approximation described 
in detail below (for the case of final results) allowed us to choose the results obtained with the 
Ax^ = 5 displacement as the most reliable. Further we adjusted the scaling factors {su) to 
improve the average quadratic approximation over the -2^ < 5 range. The necessary additional 
factors multiplying the values are shown in table 0] We see that apart from the numbers 
corresponding to the last eigenvector direction no significant adjustment is needed. 



model 


direction 


eigenvector 




"up" or "down" 


1 


2 


3 


4 


CJK 




1.04 


1.02 


1.03 


0.96 




s- 


0.96 


0.98 


0.98 


0.73 



Table 4: Multiplication factors improving the average quadratic approximation in the Zk <h 
range corresponding to all eigenvector {vi)k directions and for all {Sf} parametrizations for 
the CJK model. 

To be certain that our results are correct we need to check if the quadratic approximation 
on which the Hessian method relies is valid in the considered A^^ range. 

In the left plot of Fig. ^ we present the comparison of the dependence along each 
of four eigenvector directions (zj = 5ik for the eigenvector k) to the ideal Ax^ = zf curve. 
We see again that only the line corresponding to the eigenvector 4 does not agree with the 
theoretical prediction. Moreover it has a different shape than other lines which results from 
the scaling adjustment procedure. In the right plot of Fig. Q analogous comparison for the 5 
randomly chosen directions in the {zi] space is shown. For each of directions ^t=i -^fc ~ ^^"^ 
the ideal curve corresponds to Ax^ = z"^. In this case we observe greater deviation from the 
quadratic approximation. Finally we test the frequency distribution of Ax^ for randomly chosen 
directions in the {zi} space. Figure|21presents the frequency distribution for 1000 {zi} directions 



11 



normalized in such a way that they correspond to the ideal /S.y^ = 5, 15 or 25 {z = \/5, \/l5 
or 5 respectively). The 10, 20 and 30% deviations from the theoretical Ax^ are indicated in 
each plot. We observe that, as could be expected, the quadratic approximation worsens with 
increasing z. Still the central peak is well outlined even in the last histogram. Fractions of the 
counts within the successive lines of the increasing deviation from the theoretical predictions 
are given in table El Even for high Ax^ = 25 more than half of the counts are contained in the 
30% error range. 





10% 


20% 


30% 


Ax' = 5 


0.477 


0.652 


0.792 


Ax' = 15 


0.266 


0.498 


0.667 


Ax' = 25 


0.251 


0.401 


0.586 



Table 5: Fractions of the counts within the successive lines of the increasing deviation from the 
ideal Ax' for 1000 randomly chosen directions in the {zi} space (see text and figure 12)) for the 
CJK model. 



3.4 Collection of test CJK parametrizations 

After checking that the Hessian method gives reasonable results for the case of our CJK model 
we can create a collection of test parametrizations of the parton densities, {Sf}, for that model. 
We only need to choose the value of the magnitude t of equations (fTB|) and (fT7|). We performed 
tests with t = l,-\/5, vTO and 5 and noticed hardly any dependence of the uncertainties of 
parton densities or example physical cross-sections on its choice. We decided to apply t = 5. 
That way the {Sf} collection of eight parametrizations containing information about the CJK 
fit uncertainties was created. All sets of parameters are presented in table El 



eigenvector 


Set 




a 


P 




1 


Set 1+ 


1.916 


0.344 


0.893 


0.377 


Set 1- 


1.952 


0.258 


0.903 


0.430 


2 


Set 2+ 


1.937 


0.378 


0.841 


0.544 


Set 2- 


1.932 


0.224 


0.952 


0.271 


3 


Set 3+ 


2.462 


0.480 


0.975 


0.323 


Set 3- 


1.437 


0.129 


0.826 


0.481 


4 


Set 4+ 


1.754 


0.525 


1.912 


0.691 


Set 4- 


2.072 


0.127 


0.124 


0.185 



Table 6: Parameters of the collection of CJK parametrizations. 

Next, in table [71 we present again the parameters of the CJK model with their errors cal- 
culated within the Hessian quadratic approximation for the standard requirement of Ax' = 1. 
These uncertainties should be compared with the errors calculated by MiNUlT and shown in 
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table n They are of the same order but shghtly smaller. Note that the uncertainties given in 
table [7| can be simply multiplied by the tolerance parameter T to obtain errors which would 
correspond to an assumption of a higher Ax^ value (see Eq. [THjl . 



model 




a 






CJK 


±.JO^„0.103 


n 2qq+oo6i 


0.898^:°:?°^ 


n 404+0-066 

U.^U^_Q 054 



Table 7: The parameters of the fits for CJK model with errors calculated in the Hessian 
quadratic approximation for the standard requirement of Ax^ = 1 

All test parton distributions along with F2{x, Q"^) are further parametrized on the grid. The 
resulting FORTRAN program can be found on the web-page |77] . 

3.5 Uncertainties of the CJK parton densities 

In this section we discuss the uncertainties of the CJK model parton densities. 

In Figures OHZI the quark and gluon densities calculated in FFNScj/r models and GRV LO 
121], GRS LO [TUj and SaSlD [TT] parametrizations are compared with the CJK predictions. 
We plot for Q2 = iq and 100 GeV^ the g^(Other model)/gT(CJK) and 

^'''(Other parametrization)/g'''(CJK) ratios of the parton q'^ densities calculated in the CJK 
model and their values obtained with other models and parametrizations. Solid lines show the 
CJK fit uncertainties for Ax^ = 25 computed with the {5"^^} test parametrizations. 

First we notice that there is only one range of x, namely 0.01^a;^0.1 at = 10 GeV^ 
where the up- and down-quark densities predicted by the FFNScji<- 1 fit go slightly beyond 
the uncertainty bands. Apart from that the predictions of FFNScjk models in the case of 
all parton distributions lie between the lines of the CJK uncertainties. That indicates that 
the choice of Ax^ = 25 agrees with the differences among our three models. Moreover the 
GRV LO parametrization predictions are nearly contained within the CJK model uncertainties. 
That is not the case only for the heavy-quark densities. To hold the GRS LO and SaSlD 
parametrization curves in that range would require a much increased Ax^ value. Especially the 
SaSlD results differ very substantially from the CJK ones. 

We observe that the up-quark distribution is the one best constrained by the experimental 
data. As could be expected the greatest uncertainties are connected with the gluon densities. In 
the case of u'^ the Ax^ = 25 band widens in the small x region. Alike in the case of other quark 
uncertainties it shrinks at high-x. On contrary the gluon distributions are least constrained at 
the region of x ^ 1. That results from the fact that the gluon density contributes to F2{x, Q^) 
mainly through the 'j*G —>■ qq process which gives numerically important results only at small-x 
and therefore only in that region experimental data can constrain it. Finally we see that all 
uncertainties become slightly smaller when we go to higher from 10 to 100 GeV^. 
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4 Lagrange method for the uncertainties of the parton 
distributions 



The Hessian method, discussed above, is a very useful tool as it allows for computing of the 
parton density uncertainties in a very simple and effective way. If we find that our Ax^ 
was assumed too rigorously or on contrary too conservatively we can obtain the uncertainties 
corresponding to any other Ax^ value by simple scaling of the previous results. Still, the 
Hessian method relies on the assumption of the quadratic approximation, which in the case of 
our analysis is not perfectly preserved. Therefore, it is very important to perform a cross-check 
of the obtained results by comparing them to the corresponding results derived in a different 
statistical approach. 

There exist another method called the Lagrange multiplier method (the Lagrange method in 
short) which allows to find exact uncertainties independently of the quadratic approximation. 
This method has also been applied to the proton structure case by the CTEQ Collaboration 
jH] and the MRST group jH]. Here we utilize it to perform tests of the reliabihty of the 
Hessian approach results in the case of the CJK parton distributions. 

In the Lagrange method we make a series of fits on the quantity 

Fix, {a,}) = x\M) + AX({a,}), (19) 

each with a different but fixed value of the Lagrange multiplier A. As a result we obtain a set of 
points (x^(A),X(A)) which characterize the deviation of the physical quantity X from its best 
value Xq for a corresponding deviation of the structure function global fit from its minimum 
= X^{X) — Xo- In each of these constrained (by the A parameter) fits we find the best 
value of X and the optimal x^- For A = we return to the basic structure function fit which 
gave us {a°} parameters and allowed to calculate the best Xq value. The great advantage of 
this approach lies in the fact that we do not assume anything about the uncertainties. There 
is neither the quadratic nor any other approximation in that case. The large computer time 
consuming of the process of the whole series of minimalizations is its huge disadvantage. 



5 Examples of cross-section uncertainties in Hessian and 
Lagrange methods 

As was already said the Lagrange multiplier approach is a method on which we can rely whether 
the assumption of quadratic approximation of the Hessian is fulfilled or not. Unfortunately 
the amount of necessary computer calculations makes this approach very impractical. Still 
it can serve us as a final check of the correctness of the uncertainties that can be obtained 
using collections of our {Sf} parametrizations calculated with the Hessian method. We will 
show a comparison of the uncertainties obtained in Hessian and Lagrange multiplier methods 
calculated with the CJK parton distributions for two physical quantities. First, for the one of 
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points measured by the OPAL Collaboration [^01, as described in Section 4.3. Secondly, 
for the 7g •yq part of the 77 prompt photon production cross-section. 



We chose to exam our collection of test parametrizations first on a very simple example of Fg*^^. 
The charm-quark structure function depends only on the charm-quark and gluon distributions 
which (especially the gluon density) are not well constrained by the experimental data. On 
the other hand this dependence is not very strong in the high- a; region where the direct Bethe- 
Heitler contribution dominates. Still we expect considerable deviation of the Hessian method 
results from the correct Lagrange approach predictions. 

The jHO] OPAL analysis provided the averaged -F^c/ values in two x bins. For our purpose 
we chose the high-x point at x = 0.2 and = 20 GeV^. We calculate the CJK model 
prediction for F2^/a and its Hessian uncertainties. Further we apply the Lagrange method 
with A = ±0.33,0.66, 1, ■ ■ ■ x 10000. The results are presented in Fig. |H1 The solid line and 
crosses show the Hessian and Lagrange method predictions respectively. The dashed lines 
represent the 10 to 30% deviation from the Hessian results. As can be seen the Hessian method 
reproduces very well the Lagrange predictions in the direction of decreasing Fg^^ even for Ax^ 
greater than 100. On the other hand the agreement in the other direction corresponding to 
negative A is much worse. The lack of higher Lagrange points results from the negative A^^, 
parameter values appearing in the fits for A > 1 x 10000. 



5.2 Prompt photon production 

As the second example physical process we took a part of the prompt photon 77 production. 
Namely, for the sake of the limitation of the necessary computer time, we decided to calculate 
the direct resolved (DR) 79 — > 7g cross-section. Any cross-section of two photons can be divided 
into three parts. The direct part which is a simple 77 electromagnetic interaction, the DR being 
the part in which one photon interacts with the partons originating from the second 7 and finally 
the doubly resolved (or resolved resolved, RR) part in which both photons interact through their 
hadronic structure. Only the resolved processes contribute to the prompt photon production. 
Unfortunately calculation of the RR part is very computer time consuming. Therefore as we 
are not interested in the quantitative result of our computation but we use it only to compare 
the uncertainties obtained in two methods we omit the doubly resolved contribution to the 
cross-section. The jq 'jq process depends only on the quark distributions. We include 
the heavy-quark contributions with omitting their masses in calculations (the difference to the 
massive computation is on a 1% level). We calculate the cross-section for the photon beams 
energy = 200 GeV. That can correspond to the high energy peak of the TESLA Photon 
Collider [3T] built on the Linear Collider of the e'^e~ central mass energy of 500 GeV. 

The Lagrange method was applied with A = ±0.2, 0.4, 0.6, ■ ■ ■ x 1000. The results of the 
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comparison are presented in Fig. El Again the solid line and crosses show the Hessian and 
Lagrange method predictions respectively. The dashed lines represent the 10 to 30% deviation 
from the Hessian results. As we observe both methods agree very well in both directions of the 
change of the cross-section and in the whole range of the A^^ plotted. 

6 of the CJK fit and the data 

The CJK model presented in this article is, as described in detail in [T], a result of an improved 
analysis of j2] where our first CJKL model was given. Though we obtain a much better x^/dof 
in the case of the CJK global fit to Q^) comparing to the previous CJKL fit, its value, 

1.537, is large. Moreover, part of the improvement of the value of x^/dof is due to the exclusion 
of the TPC2^ data from the full set of the experimental F2{x,Q'^) results. Recently a similar 
global fit has been performed and presented in [20] for the case of the simple FFNS type model 
with the Bethe-Heitler process describing the heavy-quark contributions. The result of that fit 
is x^/dof = 0.938. The authors applied 134 Q^) points. In this section we try to indicate 

the main sources of the large x^/dof of our CJK fit to 182 data points. We also try to judge 
whether any exclusions of the data from the fits can be justified. 

We performed a few simple tests. First, for the case of the additional test global fit contain- 
ing the full set of the available data, mentioned in section ^TH and denoted further as FULL. 
We calculated the Xn/^n for each of the data collections as explained in 13.21 The Xn, being the 
value of the x^ of the best fit corresponding to the given experiment, divided by the number of 
the experimental points shows the agreement of each of the data collections with the global 
fit. Results are presented in the table |H1 Obviously, as different sets of the data were applied 
in this and CJK fits the corresponding numbers in tables El and |H1 are not the same. The FULL 
fit results show that the global test similarly disagree with three of the data sets, the TPC2^, 
DESY and DELPHI. 



Set 


TPC2^ 


DESY 


KEK 


ALEPH 


DELPHI 


L3 


OPAL 




2.38 


2.44 


1.20 


1.07 


2.25 


0.60 


0.97 



Table 8: The x^, being the value of the x^ of the best fit corresponding to the given experiment, 
divided by the number of the experimental points Nn for the case of the test global FULL fit 
to the full set of the available data. 

Secondly, we calculated the contribution to the x^ of the FULL fit given by each of the 208 
data points. We performed the same computation for the GRS JU] and SaSlD llj parametriza- 
tions as well as for the two FENScji^ fits presented in P (obtained without the TPC2^ data). 
We found the data which contribution to the x^ iii case of each of the mentioned fits is larger 
than 3. It appeared that there are 5 such points among 12 in the case of the CELLO exper- 
iment, 5 belonging to the DESY collection. Moreover, there are 9 such points among 22 
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in the case of the DELPHI'Ol experiment, [201; belonging to the DELPHI collection. On the 
other hand we found only 2 such points in the case of the TPC2^ data. In other experiments 
there are no more than single cases of that kind. 

As the DESY and DELPHI sets are the collections of various measurements (in the DELPHI 
experiment case they differ by the year of publication) we checked in more detail the x^/Nn 
for each of those measurements in the case of the FULL fit. The results, presented in the table 
13 confirm that as could be expected from the above test the CELLO and DELPHI'Ol results 
give main contributions to the high x^/^n of the DESY and DELPHI collections respectively. 



Set 


PLUTO 


JADE 


CELLO 


TASSO 


DELPHi'ge 


DELPHr98 


DELPHI'Ol 


Xl/Nn 


0.64 


1.12 


5.88 


1.01 


1.04 


0.10 


3.65 



Table 9: The Xn? being the value of the of the best fit corresponding to the given experiment, 
divided by the number of the experimental points for the case of the global FULL fit to the 
full set of the available data. 



Further, we performed another five test fits of the CJK model to the full data set from 
which we excluded the CELLO, DELPHI'Ol and finally pairs of CELLO, DELPHI'Ol and 
TPC2^ measurements. We denote those fits as NOCEL, NODEL (CJK fit = NOTPC) and 
NOTPCCEL, NOTPCDEL, NOCELDEL respectively. The x^ the number of experimental 
points and x^/dof for the CJK and each of the six test fits are presented in table [TUl The 
FULL fit including all available data has highest x^/dof value. The exclusion of the TPC2^ 
measurement, proposed in [2^], improves it but not very much, see the CJK fit result. The 
NOCEL and NODEL fits give much better agreement between the model and the data but 
still in their case x^/dof is larger than 1. We see that only exclusion of two of the CELLO, 
DELPHI'Ol and TPC2^ measurements at the same time leads to the x^/dof ~ 1- Finally, we 
notice that the best x^/dof is obtained in the case of the NOCELDEL fit. Central values of the 
parameters calculated in the additional test fits (not shown) and of the CJK fit (given in table 
lie inside the sum of their uncertainties. This statement is true for both cases of the CJK 
fit uncertainties: obtained from MiNOS (see table and calculated in the Hessian quadratic 
approximation (see table [7j). Only differences between the central values of the (3 parameter 
are slightly larger than the sums of the corresponding uncertainties. Therefore, we examined 
the agreement among the parton distributions computed in the CJK model in each of the fits. 
Figure ITIIl presents the comparison of the six test fit results with the CJK parton densities and 
their uncertainties obtained in the Hessian method for Ax^ = 25 and = 10 GeV^. We 
notice that the test distributions are in agreement with the CJK densities in the sense that 
they are contained within the CJK uncertainties for the assumed allowed deviation of the global 
fit from the minimum. Especially, we observe a very small difference among the CJK, FULL, 
NOCEL and NOTPCCEL parton distributions. The only deviation is observed in the case of 
the up- and down-quark densities at high-x computed for the NODEL and NOTPCDEL fits. 
The results at other values are very similar. 

We notice that the value of x^/dof for the global fits of our CJK model depends very 
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FIT 


FULL 


CJK 


NOCEL 


NODEL 


NOTPCCEL 


NOTPCDEL 


NOCELDEL 




338.7 


273.7 


262.4 


244.0 


200.7 


169.4 


179.5 


points 


208 


182 


196 


186 


170 


160 


174 


X^/dof 


1.66 


1.54 


1.37 


1.34 


1.21 


1.08 


1.06 



Table 10: The x^, the number of experimental points and x^/dof for CJK fit and six additional 
test fits. Fit denoted as FULL contains all the available data. The NOCEL, NODEL and CJK 
fits exclude the CELLO [13], DELPHI'Ol 120! and TPC2^ ^ measurements respectively. The 
NOTPCCEL, NOTPCDEL and NOCELDEL fits exclude pairs of CELLO, DELPHI'Ol and 
TPC2^ data. 



strongly on the choice of the data set applied in the fit, on the other hand, they produce similar 
parton distributions which lie within the CJK parton density uncertainties. Especially, as was 
noticed, the parton distributions predicted by our CJK and by the FULL fit, including all 
available F2{x,Q'^) data, differ very slightly. 

7 Summary 

The very first analysis of the uncertainties of the radiatively generated parton distributions 
in the real photon based on the LO DGLAP equations is given. We consider the CJK model 
presented in The estimate of the uncertainties of the CJK parton densities due to the 
experimental errors is based on the Hessian method which was recently applied in the proton 
parton structure analysis. We test the applicability of the approach by comparing uncertain- 
ties of example cross-sections calculated in the Hessian and Lagrange methods. Sets of test 
parametrizations are given, which allow for calculation of its best fit parton distributions along 
with F2{x, Q"^) and for computation of uncertainties of any physical value depending on the 
real photon parton densities. Finally, we present a detailed analysis of the of fhe CJK fit 
and its relation to the data. We show that large x^/dof of the fit is due to only a few of the 
experimental measurements. By excluding them x^/dof ~ 1 can be obtained. A FORTRAN 
program with the grid parametrization of the test parton distributions and F2{x,Q'^) can be 
obtained from the web-page P?] . 

Our work is the first trial to estimate the uncertainties of the real photon parton distri- 
butions. Future analysis of that kind should include the corrections taking into account all 
possible sources of the deviation of the experimental F2{x, Q"^) data from the photon structure 
function as described by the theoretical models. 
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z 

Figure 1: Left plot presents comparison of the x' dependence along each of four eigenvector 
directions to the ideal Ax' = zf curve. In right plot analogous comparison for 5 random 
directions in {^z^ space are shown. The ideal curve corresponds to Ax' = X] -^j' ~ -^^ 
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Figure 2: Test of the frequency distribution of Ax' for 1000 randomly chosen directions in the 
Zi space normalized so that they correspond to the ideal Ax' = 5, 15 or 25. The 10, 20 and 
30% deviations from the expected Ax' are indicated in each plot. 
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Figure 3: Up-quark densities calculated in FFNScjic models and GRV LO jSHj, GRS LO [Hlj 
and SaSlD jTTj parametrizations compared with the CJK predictions. We plot for = 10 
and 100 GeV^ the m''' (Other model) (CJK) and m''' (Other par ametrization) /«'''( CJK) ratios 
of the up-quark density calculated in the CJK model and its values obtained with other models 
and parametrizations. Solid lines show the CJK fit uncertainties for Ax^ = 25 computed with 
the set of {5*^^} test parametrizations. 
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Figure 4: The same as in Fig. IHl for the down quark. 
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Figure 5: The same as in Fig. 13 for the strange quark. 
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Figure 6: The same as in Fig. IHl for gluon. 
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Figure 7: The same as in Fig. IHl for the charm and beauty quark. Only parametrizations 
including heavy-quark distributions are compared with the CJK model predictions. 
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Figure 8: The comparison of the Lagrange and Hessian method results for the high bin F^^, 
of the OPAL 1^0] meauserement. The sohd hne and crosses show the Hessian and Lagrange 
method predictions respectively. The dashed lines represent the 10 to 30% deviation from the 
Hessian result. 
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Figure 9: The comparison of the Lagrange and Hessian method results for the direct resolved 
(DR) part of the 77 ■yq cross-section. The solid hne and crosses show the Hessian and 
Lagrange method predictions respectively. The dashed lines represent the 10 to 30% deviation 
from the Hessian result. 
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Figure 10: Comparison of the test fit parton distributions with the CJK parton densities 
and their uncertainties. Fit denoted as FULL contains all the available data. The NOCEL, 
NODEL and CJK fits exclude the CELLO [IH], DELPHI'Ol [201 and TPC2.Y |23] measurements 
respectively. The NOTPCCEL, NOTPCDEL and NOCELDEL fits exclude pairs of CELLO, 
DELPHI'Ol and TPC2^ data. We plot for = 10 GeV^ the (Other fits) /g^ (CJK) ratios of 
densities calculated in the CJK fit and their values obtained in other fits. Solid lines show the 
CJK fit uncertainties for Ax^ = 25 computed with the set of {Sf} test parametrizations. 
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